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20.  ABSTRACT 

This  report  embodies  a study  of  Ionospheric  Scintillations  as 
considered  from  the  following  aspects; 


1)  Overall  view  of  scintillations  using  a large  data  base,  and 
also  the  detailed,  individual  magnetic  storm  picture.  Solar 
flux  level,  seasonal  variation,  invariant  latitude,  and 
irregularity  region  are  discussed. 

2)  Total  Electron  Content  measurements  through  Faraday  Rotation 
and  Group  Delay  Methods.  Result_s  of  employment  of  these 
techniques  are  discussed.  The  M-factor,  for  conversion  of 
Faraday  Rotation  angles  is  studied,  and  associated  use  of  it 
is  displayed  through  a program  for  which  documentation  is 
included . 

3)  Characteristics  of  scintillations  in  the  weak  and  strong 
scatter  regimes.  Results  on  the  frequency  dependence  of 
scintillations,  power  spectra,  and  autocorrelation  of 
amplitude  fluctuations  arc  discussed. 

4)  Daytime  VHF  scintillations.  Type  1,  type  2,  and  Es^i,rregu- 
laritics  arc  discussed. 

5)  A modification  to  the  Signal  Analysis  Package,  previously 
developed,  is  described.  Application  for  use  on  the  CDC  6600 
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SCINTILLATIONS 

: 

►.  ' 

During  this  period,  work  was  done  to  focus  on  both  the 
overall  picture  of  scintillations  using  a large  data  base 
and  the  detailed,  individual  magnetic  storm  picture. 

In  conjunction  with  this  first  aim,  a high 'latitude 
model  of  scintillation  has  been  developed  using  the 
observations  of  the  existing  data  base  of  Narssarssuaq, 
Greenland,  Goose  Bay,  Labrador  and  Sagamore  Hill, 
Massachusetts.  The  available  data  consists  of: 


Narssarssuaq  9/17/68  - 9/1/74  146,700  + samples 

Goose  Bay  1/1/72  - 12/31/74  71,000  + samples 

Sagamore  Hill  12/1/69-  11/30/74  148,000  * samples 

The  model  consists  of  equations  which  yield  scintil- 
lations at  137  MHz  as  a function  of  local  time,  magnetic 
index,  solar  flux  and  month  of  the  year. 

A validation  of  this  model  was  attempted  using 
Millstone  Hill  data  obtained  by  observing  scintillations 
at  400  MHz  tracking  Transit  satellites.  The  data  base 
for  the  comparison  was  that  in  the  publication  "The 
Millstone  Hill  Radar  Propagation  Study:  Scientific  Results  " 
by  J . Evans  . 
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Since  Evans'  data  was  in  the  S4  format  at  400  Mlfz , a conversion 
using  weak  scatter  relationships  must  first  be  made  to  dB  at  137  MHz. 
This  is  plotted  in  Figure  1. 

The  measured  (400  MHz)  as  a function  of  invariant  latitude  is 

shown  in  Figure  2.  These  values,  in  three  different  Kp  groupings,  at 

the  invariant  latitude  of  Narssarssuaq,  Goose  Bay  and  Sagamore  Hill, 

along  with  the  conversion  to  dB  (137  MHz)  are  shown  below. 

NSSQ  (A  a 63°) 

S4  (400  MHz) 

1+  .0375  2.5 

3*  .05  3.6 

4.  .0775  6.3 
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Goose  Bay  (A  = 60°) 

S4  (400  MHz) 

1*  .025 

3*  .030 

4 .055 

o 

Sagamore  Hill  (A  ■ 54®) 

S4  (400  MHz) 

1 ♦ .016 

3*  .018 

4„  .021 


dB  (137MHz) 


dB  (137  MHz) 

1.5 

1.9 

4.0 


dB  (137  MHz) 
. 8 
1.0 
1.2 


This  data  used  by  Evans  was  taken  in  1972-1973  and  was  averaged  over 
local  time  to  include  all  passes.  The  model  dB  values  are  obtained 
by  using  the  local  time  equations  developed  for  each 
station  for  3 particular  Kp  values.  A solar  flux  level  of  Sf  * 108 
was  used  since  solar  flux  in  1972-73  was  of  this  magnitude.  The 

l resulting  averages  over  all  hours  and  all  months  are: 

! 

Narssarssuaq  (63°  A) 


1 

dB 

(137  MHz) 

S4  (400  MHz) 

Kp  -..5 

2.5 
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4.0 

.053 
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1.6 
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Sagamore 

Hill  (54° 
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dB 

(137  MHz) 
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.55 

.014 

' 1 

Kp  = 2.5 

.75 

.0155 

j ■ 

Kp  = 5.0 

1.23 

.022 

> 

The  model 

dB  values  (converted 

to  S^)  are 

shown  as  individual  points 

in  Figure 

2 to  facilitate  comparison  with 

Evans'  data. 
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The  results  of  Evan's  data  and  the  model  dB  values  predicted 
are  in  good  agreement  for  both  Narssarssuaq  and  Sagamore  Hill. 

' 

At  Goose  Bay,  the  model  values  are  much  lower  for  the  low 
magnetic  activity  case,  somewhat  lower  for  the  case  1(^  = 2. 5 
and  higher  for  the  case  1(^  = 5. 0. 

It  is  believed  that  the  position  of  the  trough  is  approximately 
the  invariant  latitude  of  Goose  Bay  under  quiet  magnetic 
conditions,  moving  farther  south  as  the  level  of  magnetic  activity 
increases.  The  model  data  of  Goose  Bay  reflects  this  picture. 

I 

Using  this  same  data  base,  the  seasonal  variation 
of  data  taken  at  Narssarssuaq,  Greenland,  Goose  Bay,  Labrador, 
and  Sagamore  Hill,  Massachusetts  was  also  examined.  Contour 
plots  have  been  generated  using  month  of  the  year  versus 
time  of  day  for  each  individual  station.  There  is  a further 
breakdown  of  the  data  into  low  and  high  magnetic  index  groups. 

For  the  Narssarssuaq  observations,  (Figures  3 and  4)  the 
350  km.  intersection  point  is  63“-A-  . There  arc  striking 
seasonal  contrasts  between  summer  and  winter  data  for  both 
low  and  high  Kp . The  December-February  period  shows  uniform 

1 

mean  SI  for  both  low  and  high  Kp  with  relatively  little 

diurnal  pattern  although  the  level  is  raised  for  higher  magnetic 

, activity.  The  April-June  contours  show  a strong  diurnal  pattern 

# 

with  higher  mean  SI  for  Kp=4-9.  Peak  scintillation  appears 

2300  LST  while  the  lowest  levels  are  noted  during  the  months 
of  April-June  at  08-13  LST. 


For  the  Goose  Bay  observations.  Figures  5 and  6,  the  350  km. 
intersection  point  is  60°-A-.  A higher  elevation  angle  shows 
smaller  mean  scintillation  indices  than  the  Narssarssuaq  data. 
Under  quiet  magnetic  conditions,  Kp=0-3,  very  low  mean  indices 
exist  for  the  period  of  November- January  with  little  diurnal 
variation.  Maximum  scintillation  occurs  in  the  March-May  period 
with  a definite  diurnal  pattern.  The  peak  values  occur 
^2100  LST . However,  at  Goose  Bay  the  higher  magnetic  activity 
case  is  not  simply  an  enhancement  of  the  low  Kp  case  as  occurs 

at  Narssarssuaq.  For  Kp-4-9,  minimum  SI  occurs  for  the  period 
Deccmber-February  with  a small  diurnal  variation  which  reaches 

a peak  hours  after  magnetic  midnight.  The  maximum 

activity  occurs  in  the  months  of  July-August  with  peak  values 
also  occurring  1 -2  hours  after  magnetic  midnight.  The  seasonal 
maximum  and  minimum  though  distinguishable,  are  not  as  sharply 
defined  as  at  Narssarssuaq. 


! 


The  seasonal  pattern  for  Sagamore  Hill  (at  an  intersection 
of  53°_A_  ) differs  from  that  of  the  higher  latitude  stations. 

The  irregularities  are  mainly  confined  to  the  hours  of 
darkness  for  both  low  and  high  magnetic  act ivi t y , wi th  the  peak 
daily  values  occurring  ^1-3  hours  after  magnetic  midnight , (Fig  ■ 7 £ S ) 
For  Kp=0-3,  the  seasonal  minimum  occurs  in  May-August;  for 
Kp=4-9,  it  occurs  in  April-June.  The  seasonal  maximum  appears 
to  take  place  during  the  equinoctial  months;  February- Apr i 1 
and  Scptcmber-Novcmbcr^ with  an  extension  of  the  highest 
contours  throughout  the  period  July-April  when  Kp  is  high. 
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These  Sagamore  Mil*  observations  are  felt  to  be  mid-latitude 

for  the  Kp  indices  considered  while  both  Narssarssuaq  and 

Goose  Bay  show  the  sensitivity  to  magnetic  variations  characteristic 

of  high  latitude  stations. 

For  a more  complete  view  of  scintillations,  a detailed 
picture  examining  individual  magnetic  storms  is  needed  in 
addition  to  the  overall  picture  provided  by  seasonal  contours 
and  the  model.  To  this  end,  a magnetic  storm  atlas  has  been 
proposed  including  all  the  magnetic  storms  in  the  period 
November  1 9 7 1 - Decemb er  1972,  as  noted  in  the  atlas  produced 
by  Mendillo  and  Klobuchar  (1974).  As  prototype,  data  was 
collected  for  the  magnetic  storm  of  October  18-19,  1972. 

On  October  18,  1972  at  1746  UT , a magnetic  storm  commenced. 

The  magnetic  index,  K , very  low  prior  to  the  storm,  reached 
Kp  = 3 at  SSC,  peaked  at  Kp  = 5 *v  0000  UT  October  19,  then  fell 
to  Kp=3  again  'v  1500  UT  on  October  19. 

Figure  9 is  a plot  of  the  15-minute  scintillation  index 
as  observed  at  Narssarssuaq,  Goose  Bay,  Sagamore  Hill  and 
College.  Narssarssuaq  and  Goose  Bay  exhibit  a slightly  delayed 
storm  rise  at  2030  and  2200  UT  October  18  respectively.  College 
scintillation  responds  with  medium  scintillation  ~ 0830  UT 
October  19,  while  Sagamore  Hill  shows  only  a sharp  burst  of 
scintillation  'v  2330  UT  October  18. 

In  Figure  10  is  seen  the  Narssarssuaq  and  Ottawa  magnetograms,  • 
and  the  Dst.  Both  magnetograms  show  moderate  disturbance 
levels,  peaking  in  the  period  0000-1200  UT  October  19.  The  Dst 
reaches  a maximum  value  of  -75  y at  0900  UT  October  19.  This, 


combined  with  a moderate  index,  points  to  the  irregularity 
region  developing  slowly  from  a quiet  day  picture  to  a 
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moderately  expanded  version  and  then  fading  again.  Thus,  the 
College  data  does  not  show  moderate  scintillation  activity 
until  the  irregularity  region  becomes  developed  enough  to 
reach  the  latitude  of  its  intersection  point.  Also,  the 
College  SI  seems  to  show  a seasonal  response  pattern,  with 
the  equinox  response  to  magnetic  storms  differing  from  the 
winter  dramatic  disturbance  pattern. 

The  total  electron  content  plots  of  Nars s ar s suaq , Goose 
Bay  and  Sagamore  Hill  of  October  18-19  (Fig.  11)  reveal  only 
slightly  enhanced  values  at  SSC  on  October  18.  Nighttime  TEC 
values  show  the  increased  TEC  occurring  during  many  nights  with 
magnetic  activity.  In  addition  to  somewhat  higher  values, 
spiky  areas  of  activity  are  seen  in  the  polarization  records. 

There  is  no  clear  indication  in  these  records  of  a trough 
(DMSP  pictures  are  not  available  for  this  period).  The 
scintillation  boundary,  except  for  a brief  period  between 
2330  UT  October  18  and  0030  UT  October  19,  in  the  nighttime  along 
the  75th  meridian  lies  between  Goose  Bay  and  Sagamore  Hill. 

This  short  period  is  seen  ir.  the  contour  map(Figure  12)  just 
prior  to  2400  October  19.  At  the  same  time,  the  50%  scintillation 
contour  reaches  a minimum  of  ~ 56°A  rJ  0400  UT  October  19.  Perhaps 
the  trough  lies  between  the  latitudes  of  Goose  Bay  and  Sagamore 
Hill,  but  the  synchronous  satellite  records  cannot  be  used  to 
determine  the  trough  position  in  this  case.  The  increase  in 
scintillation  index  at  Sagamore  Hill  at  2330  UT  on  October  IS  is 
accompanied  by  an  increase  in  TEC  on  the  Sagamore  Hill  TEC  records. 
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In  addition  to  the  local  time  component  of  scintillation 
shown  by  magnetic  storms,  a storm  time  component  can  be  seen 
in  the  scintillation  increases  occurring  during  some  magnetic 
storms.  Data  taken  at  various  stations  during  the  magnetic 
storms  of  August  1972  illustrate  this. 

For  large  magnetic  storms  when  the  plasmapause  moves 
well  equatorward  of  its  usual  quiettime  location,  a stormtime 
component  of  scintillations  is  recognizable  at  many  midlatitude 
stations  with  invariant  latitudes  <60°.  Figures  13  and  14 
show  the  Dst  index  and  scintillations  observed  at  stations  of 
varying  invariant  latitudes.  The  Dst  index  reached  peak  values 
''--IOOy  on  August  4-6  and  August  9.  The  high  latitude  stations 
exhibit  high  SI  values  during  these  periods  but  the  stormtime 
component  is  more  noticeable  at  the  lower  latitude  stations 
of  Sagamore  Hill,  Aberystwyth  and  Lindau.  Even  though  these 
stations  are  separated  6 hours  in  local  time,  all  show 
scintillation  peaks  % 0400  UT  on  August  4,  0000  UT  on  August  S 

and  0500  UT  on  August  6.  These  occur  at  times  of  Dst  index 
increases  and  large  TEC  depletions  at  Sagamore  Hill  (Mendillo 
and  Klobuchar,  1974).  On  August  9 M200  UT , a similar  peak 

occurs  at  these  stations  even  though  the  local  time  of  the 
stations  is  in  the  morning-noon  sector. 

The  joint  study  of  the  general  patterns  of  scintillation, 
when  combined  with  the  individual  magnetic  storm  analysis  should 
yield  a more  accurate  picture  of  scintillation  development  and 
occurrence . 
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Fig. 13  Ost  and  Scintillation  Index  at  Sagamore  Hill, 
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and  Lindau  on  3-10  August  1972. 
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II. 

Total  Electron  Content  Studies 


1 . Introduction 

The  total  electron  content  (TEC)  is  defined  as  the  total  num- 
ber of  free,  thermal  electrons  in  a vertical  column  of  unit  area  of 
ionosphere  from  the  ground  to  a height  well  above  the  peak  of  ioni- 
zation. The  use  of  artificial  satellite  beacons  permits  continuous 
monitoring  of  the  integrated  electron  content  of  the  ionosphere. 

This  is  done  by  measuring  the  effect  of  the  ionosphere  upon  radio 
signals  transmitted  from  the  satellite  and  received  at  ground  sta- 
tions. TEC  measurements  have  been  made  by  many  stations  around 
the  globe  for  many  years  by  employing  the  Faraday  rotation  tech- 
nique. Since  the  launch  of  the  ATS-6  satellite  in  1974  another 
method  of  obtaining  TEC  became  available.  The  Radio  Beacon  Exper- 
iment on  board  ATS-6  provided  the  first  reliable  differential  group 
delay  measurements  of  the  ionosphere.  The  plasmaspheric  electron 
content  is  the  content  above  the  height  at  which  the  Faraday  rota- 
tion of  a signal  propagating  through  the  ionosphere  is  negligible 
and  is  found  by  comparing  the  Faraday  and  group  delay  measurements. 
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2 . Faraday  Rotation  Method 

The  polarization  of  a plane  polarized  VHF  radio  wave  emitted 
from  a geostationary  satellite  rotates  continuously  as  it  propagates 
through  the  non- isotropic  ionosphere.  The  total  amount  of  polar- 
ization, ft,  encountered  between  the  satellite  and  the  receiving 
station  is  given  by. 


ft 


k 

f2 


B cos  0 sec  x N dh 
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where  K = 2.36x10  , f=  observing  frequency  in  Hertz,  B=  magnetic 

field  strength  in  gammas,  0=  propagation  angle,  x=  zenith  angle, 

2 

N=  local  electron  density  in  el/m  . If  the  quantity  BcosGsecx 
is  slowly  varying  with  time  it  may  be  replaced  by  its  mean  value 
along  the  propagation  path.  The  equation  of  total  polarization 
rotation  is  then  simplified  to. 


Q-  NF 


_ fhp 

where  M=  BcosOsecx  and  Np  = jq  Ndh.  The  height  hp  is  the  height  a- 
bove  which  the  remaining  amount  of  rotation  is  less  than  the  absolute 
experimental  error. 


Several  problems  occur  in  the  reduction  of  Faraday  rotation 
data.  One  is  the  problem  of  n-pi  ambiguities  in  the  data  which 
is  discussed  along  with  general  TEC  reduction  techniques  in  Cantor 
et  al  (1975)  . Another  is  the  assurance  of  absolute  rotation  val- 
ues where  the  orientation  of  the  transmitted  VHF  wave  with  respect 
to  the  satellite  spin  axis  must  be  known  as  well  as  the  reference 
polarization  of  the  polarimeter.  Absolute  calibration  of  all 
equipment  is  essential  to  valid  measurements.  Klobuchar  (1973) 
has  estimated  the  total  polarization  rotation  measurement  error 
to  be  not  better  than  - 10  degrees. 

3 . Differential  Group  Delay  Method 

An  alternative  method  of  measuring  the  electron  content  along 
the  line  of  sight  to  the  satellite  is  the  differential  group  delay 
technique.  This  method  has  been  reviewed  by  many  authors  includ- 
ing Evans  (1976)  and  Davies  (1976).  This  technique  measures  the 
difference  in  transit  times  of  two  signals  propagating  through  the 
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ionosphere  by  comparing  the  phase  of  the  modulation  envelope  be- 
tween the  carrier  and  its  sideband  at  two  widely  separated  frequen- 


i 
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cies.  The  ATS-6  satellite  carries  a mul ti - frequency  beacon  and 
is  an  ideal  source  for  this  type  of  measurement.  Group  delay 
measurements  using  the  ATS-6  beacon  have  been  recorded  at  Hamilton, 
Massachusetts  during  the  period  June,  1974  through  April,  1975. 

4 . Plasmaspheric  Electron  Content 

The  Faraday  rotation  effect  described  in  Section  II  is  strong- 
ly dependent  on  the  earth's  magnetic  field  which  is  inversely  pro- 
portional to  the  cube  of  the  distance  from  the  center  of  the  earth. 
Also,  the  electron  density  decreases  exponentially  with  altitude  a- 
bove  the  peak  density.  These  two  factors  together  indicate  that 
the  Faraday  rotation  method  yields  a measure  of  the  electron  con- 
tent of  the  lower  ionosphere,  Np,  and  is  insensitive  to  electrons 
encountered  in  the  upper  regions  of  the  sate  1 1 ite - s tat  ion  ray  path. 
The  group  delay  method,  however,  is  independent  of  the  magnetic  field 
and  gives  the  integrated  electron  content,  Ny,  along  the  entire  ray 
path.  The  plasmaspheric  content,  Np,  is  found  by  subtracting  Np 
from  N.p. 

The  ATS-6  satellite  allows  for  simultaneous  measurements  of 
Faraday  rotation  and  group  delay  along  the  same  ray  path.  A plot 
of  simultaneous  Faraday/group  delay  data  as  observed  from  Hamilton, 
Massachusetts  is  shown  in  Figure  i.  The  upper  plot  in  Figure  1 
shows  the  monthly  median  plasmaspheric  variation  for  November,  1974. 

A scatter  plot  (Figure  2)  for  the  same  period  can  be  used  for  sta- 
tistical studies  of  the  plasmaspheric  content.  Yeh  (1976)  has 
found  that  under  normal  geomagnetic  conditions  the  plasmaspheric 


content  exhibits  very  little  diurnal  variation  and  is  essentially 
a constant  value.  Initial  studies  of  plasmaspheric  contents  dur- 
ing geomagnetical ly  active  periods  have  produced  contradictory  re- 
sults (Soicher,  1976,  and  Davies  et  al,  1975).  Yeh  (1976)  and 
Davies  (1976)  have  reported  that  the  plasmasphere  becomes  deplet- 
ed during  geomagnetic  storms  and  is  steadily  replenished  by  plasma 
from  the  ionosphere  for  several  days  after  sudden  commercement . 
Further  studies  of  storm  time  effects  will  yield  rates  of  replen- 
ishment of  this  region. 

One  of  the  problems  encountered  in  deriving  plasmaspheric  con- 
tents is  that  Np  is  sometimes  of  the  order  of  0.1  N^,  thus  very 
accurate  NT  and  Np  values  are  required.  A main  factor  contribut- 
ing to  uncertainies  in  Np  is  due  to  variations  in  the  M value  con- 
version factor.  Another  consideration  is  at  what  height  the  M 
value  should  be  computed  (Yeh,  1976).  Variations  here  are  due  to 
uncentaint ies  in  the  shape  of  the  ionspheric  Ne(h)  profile. 


5 . M-  Value  Determination 

The  M factor  is  used  to  convert  Faraday  rotation  angles  in  in- 
tegrated electron  contents.  For  an  ideal  geostationary  satellite 
in  a perfectly  synchronous  orbit  with  orbital  inclination  and  ec- 
centricity exactly  equal  to  zero  the  M value  is  indeed  a constant 
and  needs  to  be  determined  only  once.  In  real  life,  even  if  the 
ideal  orbit  is  initially  achieved  it  will  soon  decay  due  to  the 
perturbing  forces  if  the  sun,  moon,  and  variations  in  the  earth's 
gravitational  field. 

An  inclination  in  the  satellite's  orbit  results  in  a diurnal 


32 


oscillation  in  latitude  of  the  sub - s at e 1 1 i te  point  around  a mean 
value  of  latitude.  Any  ellipticity  in  the  orbit  will  produce  a 
similar  oscillation  in  longitude.  Several  examples  of  this  are 
shown  in  Figure  3 and  Figure  5.  Figure  3 is  a plot  of  the  actual 
hourly  variation  in  latitude  and  longitude  of  the  sub - s at e 1 1 i t e pos- 
ition of  I2F3  which  has  an  inclination  of  6.6  degrees  and  eccentri;- 
ity  of  .002  degrees.  The  corresponding  M values  for  an  observer 
at  Lunping,  Taiwan  are  shown  in  Figure  4.  A plot  of  the  sub-sat- 
ellite position  of  ATS-3  which  has  an  inclination  of  5.9  degrees  and 
eccentricity  of  .002  degrees  appears  in  Figure  5 and  the  correspond- 
ing M values  for  an  observer  at  Hamilton,  Massachusetts  are  shown 
in  Figure  6.  The  amplitude  of  the  diurnal  variation  in  M value 
for  1 2 F 3 as  observed  from  Lunping  is  8.8%  from  the  constant 
mean  value  which  would  traditionally  be  used  in  the  Faraday  computa- 
tion. The  amplitude  of  the  diurnal  variation  in  M value  for  ATS-3 
as  observed  from  Hamilton  is  only  1,5.5%.  There  are  several  effects 

to  be  noted:  1.)  For  a geostationary  satellite  observed 

from  a station  at  middle  or  high  latitudes 
the  relative  change  in  the  M value  with  time 
is  smaller  than  other  experimental  errors  and 
can  be  ignored.  For  observes  in  or  near  equa- 
torial regions  the  relative  change  in  M value  is 
large  and  use  of  the  traditional  constant  M value 
conversion  factor  is  invalid.  An  actual  value  cf 
M appropriate  to  the  time  of  observation  must 
be  calculated  if  one  is  to  obtain  a reasonable 

measure  of  N... 

F 

2.)  As  the  inclination  of  the  orbital  plane  of 
the  satellite  increases  with  time  the  sub-satell- 
ite position  will  have  greater  latitudinal  os- 
cillations and  thus  a larger  fluctuation  in  M 
value  will  occur. 
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A computer  program  has  been  developed  which  will  calculate 
a satellite's  position  as  a function  of  time  for  a given  set  of 
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orbital  elements  and  will  also  calculate  the  corresponding  M value 
for  an  observer  at  a given  station  as  shown  in  Table  1.  This  pro- 
gram is  described  in  more  detail  in  Section  VI.  The  program  pro- 
vides a relatively  simple  corrective  procedure  which  can  be  used 
by  the  monitoring  stations  for  obtaining  more  accurate  TEC  values. 

This  is  accomplished  by  use  of  a look-up  table,  an  example  of 
which  is  given  in  Table  2 and  Table  3,  which  converts  Faraday  ro- 
tation angles  expressed  in  the  polarimeter  polarization  channel 
strip  chart  units  of  pi's  and  millimeters  into  corrected  TEC  values. 

These  tables  have  been  designed  to  minimize  the  amount  of  time  and 
effort  that  is  required  of  an  observer.  The  conversion  tables 
are  currently  in  use  at  observing  stations  in  Athens,  Greece/  Goose 
Bay,  Labrador  / Hami 1 ton , Massachusetts  / Kennedy  Space  Flight  Center, 

Florida/  Ramey  Air  Force  Base,  Puerto  Rico/Osan,  Korea  and  Lunping^ 

Taiwan . 

Program  Support 

A satellite  ephemeris  program,  LOKANGL,  was  developed  by  the 
Analysis  and  Simulation  Branch  of  the  Air  Force  Geophysics  Laboratory. 

Under  this  Contract  (F19628-76-C-0003)  the  program  has  been 
modified  and  expanded  to  aid  in  improving  the  accuracy  of  reported 
TEC  values. 
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OUTPUT : For  any  given  satel 1 ite- station  pair  the  program  will  cal- 
culate the  satellite's  orbit  and  will  produce  the  sub- s ate  1 1 i t e and 
sub-ionospheric  latitudes  and  longitudes,  the  satellite's  right 
ascension  and  declination,  elevation,  azimuth,  and  range  as  a 
function  of  time  as  observed  from  the  user  specified  station. 

Other  output  options  include: 

(1)  Position  and  velocity  vectors 

(2)  Sub-satellite  data 

(3)  Osculating  Elements 

(4)  Table  of  M-values  as  a function  of  time 

(5)  Conversion  table  (TEC  (constant  M)  to  TEC  (varying  M)  ) 

(6)  Faraday  rotation  to  TEC  conversion  table 

Any  of  these  output  options  may  be  generated  from  one  of  the 
following  input  sets: 

(1)  A 2-card  orbital  element  data  set 

(2)  A 2-card  input  set  of  position  and  velocity  vectors 

(3)  A 5-card  SPADATS  orbital  element  data  set 

The  5-card  SPADATS  orbital  element  set  has  been  used  exclusive- 
ly and  only  this  input  set  will  be  Described  in  this  report. 

\ 

FUNCTIONAL  DESCRIPTION 

I . Ephemcris  program 

NAME 
LOKANGL 

CONVRT 
DELEM 
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NAME 

DREV 


NUMINT 

INTGND 

DENSEL 

OSCTMN 

PVTELM 

MNTOSC 

ELMTPV 

SHPRDC 

KEPECN 

SPSETC 

SPRFCO 

NEW 

CRMXA 

SPTRAN 

SPROU 

MULTIP 

JDATE 

CDATE 

WRSTA 

WRTAP 
ANGLE  1 
ANGLE  2 


1 

! 


DESCRIPTION 

obtains  change  per  revolution  in  orbital  elements 

computes  changes  per  rev  due  to  drag  by  Simpsons  rule 
integration 

supplies  NUMINT  with  integrands  for  each  integ.  step 
density  model 

converts  osculating  pos ition- velocity  to  mean  elements 
converts  pos i t ion - ve 1 oci ty  to  orbital  elements 
converts  elements  to  osculating  pos i t ion- ve 1 oci t y 
converts  elements  to  position-velocity 
computes  short  periodic  terms 
converts  mean  anomaly  to  eccentric  anomaly 

sets  constants  and  program  initializations  for  all  routines 
computes  refraction  correction  due  to  atmosphere 
converts  units  of  output  data 
prints  heading  on  each  output  page 

transformation  matrix  and  station  coordinates  routine 
prints  out  BDC  and/or  binary  ephemeris  tapes 
performs  matrix  multiplication 

converts  Gregorian  date  to  Julian  date-limits  Jan  1,1950 
to  Dec  31,1999 

converts  Julian  date  to  Gregorian  date-sa^  limits 

reads  M-value  and  TEC  conversion  table  output  opt i ons 
and  prints  ephemeris  predictions 

extracts  lookangle  predictions  from  binary  tape 

reduces  an  angle  to  interval  -ir  to  +tt 

reduces  an  angle  to  interval  0 to  2tt 
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M-VALUE  ROUTINES 


i NAME 


mgEos 


FIELDG 

FIELD 


HOUR 


HEADER 


HEAD 


DESCRIPTION 

CALLED  FROM  WRSTA,  calculates  M-value  as  function 
of  satellite  position,  calls  necessary  routines 

sets  initializations  and  performs  geometric  calcula- 
tions 

calculate  geomagnetic  field,  updates  the_POGO  coef- 
ficients to  time  of  observation,  prints  M-value  table 

calculates  hour  angle 

prints  AFGL  header  on  every  page  of  output 
prints  heading  on  M-value  table 


III.  TEC  CONVERSION  TABLE  ROUTINE 


NAME 


TABLE 


DESCRIPTION 

produces  the  following  output 

(1)  TEC  (constant  M)  to  TEC  (varying  M) 
conversion  table 

(2)  Faraday  rotation  to  TEC  conversion  table 


INPUT 


The  input  which  is  required  by  program  LOKANGL  is  of  two  forms 
One  is  a series  of  input  cards  which  are  to  be  submitted  with  the 
program  deck.  The  other  is  an  input  tape  (tape  2)  on  which  the 
POGO  (8/69)  geomagnetic  field  coefficients  reside.  A sample  of 
the  control  deck  is  included  later  in  this  report. 
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I.  INPUT  CARD  FORMAT 
CARD 


NO. 

COLUMN  FORMAT  UNITS 

1 

1 

2 

1 

1 

4-7 

A4 

46-48 

12 

49-51 

F 3 . 0 

i 

3 

1 

9-13 

F5 . 0 

23-31 

F9 .4 

degrees 

32-40 

F9 .4 

degrees 

41-49 

F9 . 4 

degrees 

50-58 

F9 . 4 

59-67 

F 9 . 4 

degrees 

4 

9-19 

Ell  . 8 

revs/ day 

21-31 

El  1 . 9 

rev/day 

33-40 

E 8 .5 

deg/day 

42-49 

E 8 . 5 

deg/day 

51-58 

E8 . 5 

60-67 

E 8 . 5 

deg/day 

5 

: 

1 

l 

9-19 

El  1 . 8 

rev/day'* 

21-31 

Ell  .8 

rev/day4 

33-40 

E8 .5 

2 

rev/day 

42-49 

E8.5 

rev/day^ 

DESCRIPTION 

"5 "=  5 card  element  set  to  be  used 

"1"=  first  card  of  element  set 

satellite  name 

year  of  orbital  elements 

Julian  day  number  of  orbital  elements 

"2"=  second  card  of  element  set 

Epoch=  fractional  part  of  orb  element 
date'=0.0  if  epoch=  Ohr  U.T. 

Mean  anomaly 

Right  Ascension  of  ascending  node(H) 
Argument  of  Perigee  (w) 

Eccentricity  (e) 

Inclination  (i) 

Mean  motion^  1440.  (mins/day) /period 
(mins/rev)=  n 

n/2 

h 

0) 

e/day 

« 

i 

"4"  = fourth  card  of  orbital  element 
n/6 
■h’/2A 

h‘/2 

w/2 


set 
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CARD 

NO. 

COLUMN 

FORMAT 

UNITS 

DESCRIPTION 

6 

1 

"5"  = fifth  card  of  orbital  ele- 
ment set 

9-19 

Ell  .8 

E . R . 

Semi-major  axis  (a)  E.R.=earth 
radi  i 

1 l 

1 

21-30 

E10 . 7 

E . R . /day 

a 

32-41 

E10 . 7 

E . R . / day ^ 

a/ 2 

f 7 

1-3 

13 

number  of  stations  requesting  look 
angle  output 

4-6 

13 

if  number  of  stations  1,  code  in 
column  6:  0=  output  as  function  of 
time 

1=  output  as  function  of 
time 

8 

station  parameters 

1-5 

15 

I . D . of  station 

8 

13 

"2"  = geodetic  system 

9-23 

E 1 5 . 8 

degrees 

station  longitude  ( + W) 

24-38 

E 1 5 . 8 

degrees 

station  geodetic  latitude  (+N) 

39-53 

El  5 . 8 

meters 

station  height 

61-72 

2A6 

name  of  station 

9 

start  and  stop  times  of  ephemeris 
table 

1-15 

E 1 5 . 8 

seconds 

print  interval  (t) 

L 

17-18 

12 

start  month 

20-21 

12 

day 

! 

23-24 

12 

year  (last  two  digits) 

26-29 

F4  . 1 

hour 

V 

31-34 

F4 . 1 

minutes 

36-41 

F 7 . 3 

seconds 

43-44 

12 

end  month 

46-47 

12 

39 

day 

L — 

— ^ 

t t- 

-m--  ^ 

CARD 

NO. 

COLUMN 

FORMAT 

UNITS 

DESCRIPTION 

9 

49-50 

F 4 . 1 

end  year  (last  two  digits) 

52-55 

F 4 . 1 

hour 

57-60 

F 4 . 1 

minutes 

62-67 

F 4 . 1 

seconds 

10 

print  codes  all=0  when  look  angles 
requested 

2 

12 

=0;  =1  to  print  position/velocity 
vectors 

4 

12 

= 0;  =1  to  print  s ub - s ate  1 1 i t e position 

6 

12 

=0;  =1  to  print  osculating  elements 

8 

12 

=0;  =1  to  print  station  observation 
data 

9-10 

12 

=0;  =1  to  write  binary  tape 

11-13 

E 1 5 . 8 

=1.0;  single  correction  of  range  for 
refraction 

11 

16 

11 

output  options  for  M-value  and  Con- 
version Tables 
"0"=  M-value  table  only 
"1"=  M-value  and  Faraday  to  TEC  Con- 
version Tables 
" 2"=  all  tables 

33-38 

F6 . 2 

MHz 

operating  frequency 

5 2-55 

F4  . 1 

MM 

Polarization  reference  va 1 ue  (mi  1 1 i me t e r s 

60-63 

14 

Year  of  requested  M-value  prediction 

1 2 

1-3 

"END" 

5-6 

"OF" 

8-14 

"PROBLEM" 

13 

1 

"9"  = end  of  all  input 

II.  INPUT  TAPE  FORMAT 

The  POGO  (8/69)  set  of  geomagnetic  field  coefficients  must  be 
attached  to  the  local  file  TAPE2  before  program  execution.  The  for- 
mat for  this  file  is  as  follows; 


CARD  NO. 

COLUMN 

VARIABLE  NAME 

FORMAT 

DESCRIPTION 

1 

1 

J 

11 

=0  if  coefficients  were  de- 
rived using  oblate  earth 
/0  if  coefficients  were  de- 
rived using  spherical  earth 

1 

2 

K 

11 

=0  if  coefficients  are  Schmidt 
normalized 

/0  if  coefficients  are  Gauss 
normalized 

1 

4-9 

TZERO 

F6 . 1 

epoch  for  which  coefficients 
were  derived  = 1960.0 

1 

10-72 

AID 

1 0A6 jA3 

Alphabetic  ident i f i ca t i cn 

Each  succeeding  card  contains  the  variables  N , M , GNM , HNM,  GTNM , 
HTNM,  GTTNM,  HTTNM,  in  the  following  format  (2I3.6F11.4)  where  GNM, 
HNM,  GTNM,  HTNM,  GTTNM,  and  HTTNM  are  the  spherical  harmonic  coeffi- 
cients in  units  of  gamma  (10  ^ gauss). 

x,  , „ - m 


N = n + 1 

GTNM  = 

gm 
n 1 

M = m+ 1 

HTNM  = 

1>  , 

^ time  derivatives 

n 

GNM  = gm  ) 
" \ 

spherical 

GTTNM 

E C 
* 

ii 

at  the  epoch  da:e 

I 

HNM  = hm 
n 

harmonic 

coefficients 

HTTNM 

VI 
= h 

n 

tUUllU,iCHLD  II 


Termination  of  the  coefficient  set  is  accomplished  by  a flag 
card  where  N is  set  equal  to  zero.  A list  of  the  POGO  (8/69)  coeff- 
icient file  is  given  in  Table  4 of  this  report. 
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EXECUTION  REQUIREMENTS 

Core  memory  required  = 110  K octal 

Compile  time  = 17.9  decimal  seconds 

Execution  time  = 8.7  seconds  per  station/month 

Since  this  program  will  be  used  to  produce  output  ta- 
bles for  many  stations  each  month  one  would  like  to  minimize 
the  execution  requirements.  One  way  of  accomplishing  this 
is  to  store  the  binary  object  code  as  a us er- 1 ibrary , thus 
avoiding  the  need  for  compiling  the  program  each  time  it  is 
use  while  greatly  reducing  the  operating  costs. 


TAPE  REQUIREMENTS 


The  following  tapes  must  be  included  on  the  program  card: 


TAPE  1 

output 

tape 

containing 

the  M-valuc  and  TEC  tables 

TAPE  2 

input  tape  containing 

the  P0G0 

(8/69}  coefficient 

TAPE3 

i 

scratch 

t ape 

, TAPE4 

scratch 

tape 

TAPE5 

input  ; 

card 

reader 

TAPE6  : 

■ 

output ; 

line 

printer  ; 

contains 

ephemeris  data 

* 
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USE  OF  PROGRAM 


In  order  to  use  the  Faraday  to  TEC  conversion  correctly  an 
observer  should  first  verify  that  the  station  location,  operating 
frequency,  satellite  being  observed,  and  reference  polarization 
value  match  those  listed  on  the  computer  table  headings. 

The  observer  then  turns  to  the  page  of  the  conversion  tables 
with  the  proper  pi-value  and  finds  the  millimeter  reading  which 
matches  the  polarization  channel  chart  reading.  He  then  scans  across 
the  table  to  find  the  final  TEC  value  in  the  column  which  corresponds 
to  the  time  of  the  observation. 

For  example,  referring  to  the  sample  output  in  Table  3 one  can 

easily  see  that  for  an  observer  at  Lunping,  Taiwan  monitoring  the 

satellite  IS2F3  which  has  a nominal  longitude  of  -161.9  degrees  east, 

on  a frequency  of  136.44  MHz,  with  a reference  polarization  value  of 

13.5  millimeters  during  the  month  of  September,  1976,  a polarization 

chart  reading  of  2n,  15  millimeters  at  14  hrs  universal  time  corre- 

16  2 

sponds  to  a TEC  value  of  13.8  x 10  el/m  . 
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SAMPLE  CONTROL  DECK 


0JL0K/T90/C1-U  lOOOO. 

ATTACH  / JL IB i JMJMBX2  93  9/  ID=  J0HAM50N. 
ATTACH/ TAPK?. /TOG 0X2 93  9/  ID-  JOHANS ON. 

l ion  Aim  jlid ) 

LDSET (PRESET  = ZERO ) 
l.OKAIvGL. 

P.EV'  IKD/  TAPE  1 . 

C OP Y, TAPE  1 /OUTPUT. 

EXIT. 

7/8/9 

5 

1 12  F3 

2 0.0  55.8347  54. 

.3  1 .002989  / -.01 

» 

4 

5 6.6096106 

II 

47  .-121.17  25.0 

3600.  09  15  76  00.0  00.0  0 

0 0 0 0 0 1 . 0 

OUTPUT  OPT  1 CM  = 1 FREOUENC Y =13  6. 

1T.KD  OF  PROBLEM 


2535 


JOHANS ON 


-121.17 


076127 

55.8347  54.0444  298.4900  .0019999 
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SIGNAL  STATISTICS  0?  IONOSPHERIC  SCINTILLATION 
IN  THE  WEAK  AND  STRONG  SCATTER  REGIMES 

1.  INTRODUCTION 

For  a successful  design  of  ground  to  space  communication 
links,  a knowledge  of  ionospheric  scintillations  at  a desired 
frequency  as  well  as  the  fine  structure  of  the  fluctuations 
is  required.  At  mid- lat itudes  where  the  amplitude  fluctuations 
are  small,  ground  based  scintillation  measurements  in  conjunction 
with  the  we  1 1 - de ve 1 oped  weak  scattering  theories  allow  predictive 
systems  of  scintillations  to  be  successfully  built.  On  the  other 
hand,  in  the  equatorial  region,  where  strong  irregularities  extend 
over  several  hundred  kilometers,  amplitude  fluctuations  as  large 
as  25  dB  in  the  VHF  range  and  even  2 dB  at  frequencies  as  high  as 
6 GHz  are  often  observed.  The  deleterious  effects  of  such  intense 
scintillations  on  communication  links  are  indeed  very  formidable. 

In  view  of  the  paucity  of  experimental  data  defining  the  structure 
of  ionospheric  scintillations  in  the  strong  scatter  regime  and  the 
fact  that  scattering  theories  in  this  domain  are  not  fully  developed, 
it  has  not  been  possible  to  develop  a successful  predictive  system 
of  scintillations  in  the  equatorial  region  for  use  in  commun ic  at  ion 
engineering. 

In  this  report,  we  present  such  a study  of  the  characteristics 
of  ionospheric  scintillation  in  the  weak  and  strong  scatter  regimes. 
Results  on  the  frequency  dependence  of  scintillations,  power 
spectra  and  autocorrelation  intervals  of  amplitude  fluctuations  in 
the  two  scattering  domains  are  discussed. 
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OBSERVATIONAL  MATERIAL  AND  METHOD  OF  ANALYSIS 


2 . 

The  study  has  been  based  on  simultaneous  amplitude  scin- 
tillation measurements  of  137  and  360  MHz  transmissions  from 
ATS-6  carried  out  at  Huancayo,  Peru.  The  data  used  in  the 
study  correspond  to  a high  level  of  scintillation  activity 
and  were  recorded  over  the  interval  2215-2400  LT  on  Nov.  1,  1974 
and  2018-2333  LT  on  Nov.  4,  1974.  In  order  to  study  the  behavior 
of  the  wavelength  dependence  of  scintillation,  low  fluctuation 
level  scintillation  data  obtained  at  Sagamore  Hill  Radio 
Observatory,  Hamilton,  Mass.,  have  also  been  used.  Sagamore 
Hill  data  were  obtained  from  ATS-5  transmissions  at  137  and 
412  MHz  on  April  12  and  May  18,  1973.  The  data  were  recorded 
on  FM  analog  tape  which  was  later  digitized  at  the  rate  of 
6 s amp  1 e s / se cond . From  each  15-minute  data  block  from  Huancayo 
and  every  5-minute  block  from  Sagamore  Hill^the  index  of 
scintillation  has  been  derived  which  is  defined  as  the  second 
moment  of  intensity  normalized  to  the  average  intensity. 

The  137  and  360  MHz  scintillation  data  recorded  at  Huancayo 
have  also  been  processed  to  obtain  the  autocorrelation  intervals 
and  the  power  spectrum  of  scintillations  by  the  use  of  FFT 
al gorithm . 


56 


3.  FREQUENCY  DEPENDENCE  OF  SCINTILLATION 
IN  THE  WEAK  AND  STRONG  SCATTER  REGIMES 

From  channel  modelling  considerations  there  is  considerable 
interest  in  the  study  of  frequency  dependence  of  scintillation 
as  it  provides  an  estimate  of  scintillation  at  a desired  frequency 
by  extrapolating  scintillation  measurements  at  another  frequency. 
Recent  in-situ  measurements  at  F-region  heights  indicate  that 
the  power  spectral  density  of  the  irregularity  wavenumber  spectrum 
exhibits  a power  law  form  with  one-dimensional  spectral  index 

2.  The  dependence  of  three-dimensional  irregularity  power 
spectral  density  $.,(10  on  the  irregularity  wavenumber  k may  thus 


be  represented  as; 


*N(k)  = k-p 


where  p = + 2 n,  4 for  F-region  irregularities. 

It  has  been  shown  that  in  the  weak  scatter  limit  or  small 
amplitude  fluctuations,  the  above  form  of  irregularity  spectrum 
yields  a frequency  dependence  of  index  of  scintillation  as; 

S4  f'(P  + 2)/4  <r  f‘n  (2) 

Thus  the  exponent  n for  the  frequency  dependence  of  scintillation 
is  related  to  the  spectral  index  p of  three-dimensional  irreg- 
ularity wavenumber  spectrum  as; 


For  F-region  i rregul a ri t i es  if  we  assume  that  pRs4  the  exponent 
for  wavelength  dependence  q is  expected  to  have  a value  of  about 
1.5  for  weak  scintillations. 

Based  on  Sagamore  Hill  observations  of  137  and  412  MHz 
scintillations  of  ATS-5  transmissions  on  April  12  and  May  18, 
1973,  we  obtained  the  exponent  n.  In  order  to  examine  if  the 
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weak  scatter  conditions  are  violated  we  studied  the  variation 
of  n with  index  at  137  MHz,  corresponding  to  the  lower  limit 

of  the  frequency  pair.  Most  of  the  observation  points  pertained 


to  < 0.4  at  137  MHz  conforming  to  weak  scatter  limit.  The 

. 

j data  of  April  12  , 1973  yielded  a median  value  of  i)  ^ 1.2  while 

May  18,  1973  data  yields  q 1.5.  These  values  of  q correspond 
to  spectral  indices  varying  between  3.8  to  4.0  which  is  very 
close  to  that  obtained  from  in-situ  measurements  of  F-region 
irregularity  power  spectrum. 

The  frequency  dependence  of  scintillation  in  the  regime  of 
strong  amplitude  fluctuations  was  studied  by  the  use  of  scin- 
tillation measurements  at  Huancayo  with  ATS-6  transmissions  at 
137  and  360  MHz.  It  may  be  observed  that  the  exponent  n 
progressively  decreases  with  increasing  S4  in  the  strong 
scatter  regime.  For  > 0.9,  q % 0 , indicating  thereby  that 
scintillations  become  invariant  of  wavelength  for  very  strong 
scattering. 

Numerical  solutions  to  the  problem  of  multiple  scattering 
indicate  that  scintillations  gradually  saturate  as  scattering 
conditions  become  stronger.  In  the  case  of  our  two  frequency 
(137/360  MHz)  observations,  scintillations  saturate  first  at 
lower  frequency  while  the  scintillation  index  of  the  higher 
frequency  still  registers  3n  increase  with  increasing  strength 

> 

of  scattering.  This  explains  the  observed  progressive  decrease 
of  q index  with  S^.  Finally,  when  the  scattering  becomes  so 
strong  as  to  drive  both  137  and  360  MHz  to  the  saturated  scin- 
tillation regime,  frequency-independent  scintillation  (q  0) 
should  result.  The  experimental  results  may  thus  be  explained 
in  the  framework  of  multiple  scatter  theories. 
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4.  POWER  SPECTRUM  SCINTILLATIONS  IN  THE  WEAK 
AND  MULTIPLE  SCATTER  REGIMES 

The  power  spectrum  of  scintillations  in  the  equatorial 
region  was  studied  and  its  variation  with  scintillation  activity 
was  examined. 

Figures  2 and  3 illustrate  the  charac t eri s t i cs  of  power 

spectrum  of  scintillations  at  137  and  360  MHz  over  two  15-minute 

intervals.  The  spectra  in  Figure  2 were  obtained  when  the  S4 

indices  at  137  and  360  MHz  were  0.36  and  0.13  respectively  and 

the  case  therefore  corresponds  to  weak  scatter  conditions. 

The  noise  and  frequency  resolution  error  bars  as  computed  at 

0.1  Hz  are  indicated  in  the  diagram.  The  spectral  structure 

around  1 Hz  corresponds  to  satellite  spin  modulation.  It  may 

be  observed  that  the  spectra  are  characterized  by  a relatively 

flat  portion  and  then  a linear  roll-off  portion.  The  spectral 

slopes  of  the  roll-off  portions  correspond  to  f 3.6-0. 5 an{j 
- 2 3±0  5 

f . If  we  consider  the  spectral  slope  to  be  approximately 

3,  then  from  weak  scatter  theory  an  irregularity  power  spectral 
index  of  4 is  obtained  in  agreement  with  in-situ  measurements 
of  F region  irregularities 

Figure  3 illustrates  the  power  spectrum  of  scintillations 
when  multiple  scatter  conditions  prevailed  at  both  frequencies, 
the  indices  at  137  and  360  MHz  being  0.95  and  0.84  respec- 
tively. The  marked  variation  of  the  spectral  features  in  this 
case  as  compared  to  tha$  under  weak  scatter  domain,  illustrated 
in  the  previous  diagram,  is  evident.  The  flat  portion  of  the 
spectrum  extends  to  much  higher  frequency  and  the  slopes  of  the 
roll-off  portions  become  very  steep  as  strong  scattering  is 
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encountered.  The  frequency  extension  of  the  flat  portion 
indicates  that  the  noise  band  becomes  more  white  in  the  strong 
scatter  case  and  is  more  red  in  the  weak  scatter  condition. 

A close  examination  reveals  that  the  limiting  frequency 
of  the  flat  portion  is  higher  at  137  MHz  as  compared  to  360  Mllz. 

The  results  were  confirmed  independently  by  determining  directly 
the  autocorrelation  intervals  at  the  two  frequencies.  That  the 
autocorrelation  interval  is  shorter  at  lower  frequency  as  compared 
to  the  higher  frequency  is  opposite  to  what  is  expected  under 
weak  scatter  theory. 

From  a knowledge  of  irregularity  velocity,  the  autocorrelation 
intervals  may  be  transferred  to  autocorrelation  distance.  This 
result  has  important  implications  in  space  diversity  systems. 

We  find  that  for  combating  weak  scintillations  the  space  diversity 
at  lower  frequencies  is  to  be  larger  than  that  required  at  higher 
frequency  but  the  reverse  situation  will  hold  when  strong  scin- 
tillations are  encountered. 


60 


LIST  OF  FIGURES 


Figure  1 


Figure  2 


Figure  3 


The  variation  of  the  power  law  exponent,  n, 
defining  the  frequency  dependence  of  scintillations 
(f  ^)  , with  S^  index  at  137  MHz. 

The  power  spectra  of  amplitude  scintillations 
at  137  and  360  MHz  when  scattering  is  weak. 

The  power  spectra  of  scintillations  at  137  and 
360  MHz  when  scattering  is  strong. 
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ON  THE  NATURE  OF  ELECTROJET  IRREGULARITIES 
RESPONSIBLE  FOR  DAYTIME  VHF  SCINTILLATIONS 

1.  INTRODUCTION 

Scintillation  measurements  made  at  Huancayo,  Peru  by  the  use 
of  41  and  140  MHz  transmissions  of  ATS-6  indicate  that  scintilla- 
tion is  almost  a regularly  occurring  phenomenon  during  daytime 
near  the  magnetic  equator.  In  view  of  the  much  lower  level  of 
daytime  scintillation  as  compared  to  nighttime,  the  character- 
istics of  daytime  scintillation  have  remained  unexplored  whereas 
equatorial  nighttime  scintillation  has  been  extensively  studied. 

The  most  viable  portion  of  the  equatorial  ionosphere  where  plasma 
density  irregularities  are  generated  during  daytime  corresponds 
to  the  electrojet  region  which  is  to  be  associated  with  the 
observed  daytime  scinti 1 lations . 

The  nature  of  electron  density  irregularities  with  meter 
wavelength  present  in  the  equatorial  electro  jet  has  been  exten- 
sively explored,  experimentally,  by  the  radar  spectral  techniques 
and  theoretically  in  the  light  of  plasma  instability  mechanisms. 
These  studies  reveal  the  existence  of  two  distinct  types  of  meter 
wave  irregularities,  namely,  types  1 and  2.  Type  1 irregularities 
are  considered  to  be  generated  by  a modified  two-stream  instability 
when  the  electron  drift  velocity  relative  to  the  ions  exceeds  the 
ion-acoustic  velocity  of  about  360  m/sec.  Type  2 irregularities, 
on  the  other  hand,  are  generated  by  gradient  drift  instability 
mechanism  when  the  electron  drift  velocity  exceeds  a velocity  of 
about  30  m/sec  in  a region  of  steep  density  gradient. 

In  view  of  the  fact  that  a backscatter  radar  is  sensitive  to 
half  the  transmitted  wavelength,  sweep  frequency  radars  have  been 
operated  over  the  frequency  range  of  IS  ISO  MHz  to  study  the 


electrojet  irregularities  with  wavelengths  ranging  between  10m  - 
Ira.  These  studies  indicate  that  the  long  wave  irregularities 
correspond  to  type  2 generated  by  gradient  drift  instability 
mechanism  whereas  short  wave  irregularities  correspond  to  type  1 
being  generated  by  two-stream  instability  mechanism.  In  an 
effort  to  extend  the  study  of  electrojet  irregularities  to 
wavelengths  longer  than  those  explored  by  radar  measurements, 
the  features  of  equatorial  sporadic-E  (E  ) in  relation  to  radar 
measurements  have  been  recently  explored.  The  sporadic-E  signature 
on  ionograms  corresponds  typically  to  30-80m  irregularity  wave- 
lengths. It  has  been  found  that  the  long  wave  Eg^  irregularities 
are  related  to  type  2 meter  wave  irregularities  which  are 
generated  by  the  gradient  drift  mechanism.  This  indeed  conforms 
to  the  theoretical  predictions.  Recent  theoretical  studies  show 
that  the  gradient  drift  mechanism  generates  still  longer  irreg- 
ularity wavelengths  in  the  kilometer  range.  From  experimental  as 
well  as  theoretical  standpoint  there  exists  considerable  interest 
in  exploring  the  long  wavelength  end  of  electrojet  irregularities. 

The  present  report  based  on  VIIF  scintillation  measurements 
and  ionograms  obtained  at  Fuancayo  as  well  as  NOAA  oblique  VHF 
radar  measurements  made  at  Jicamarca  attempts  such  a study  of  the 
characteristics  of  electrojet  irregularities.  It  is  known  that 
scintillations  arise  from  irregularities  with  wavelengths  of  the 
order  of  Fresnal  dimension  of  the  exploring  wavelength  at  the 
irregularity  height.  For  ATS-6  observations  from  lluancayo  the 
Fresnel  dimensions  at  41  and  140  MHz  with  110km  slant  range  of 
the  electrojct  region  correspond  to  1km  and  0.5km  respectively. 

The  study  of  E on  ionograms  obtained  at  lluancayo  provided 
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simultaneous  information  on  30-80m  irregularity  wavelengths  whereas 
50  MHz  oblique  radar  backscatter  measurements  carried  out  at 
Jicamarca  reveal  the  existing  types  of  irregularities  at  3m 
wavelength.  This  study,  therefore,  covers  the  wide  range  of 
irregularity  wavelength  from  a few  meters  to  about  a kilometer  in 
the  equatorial  electrojet. 

2.  RESULTS  AND  DISCUSSIONS 

On  normal  days,  the  onset  of  scintillations  occur  immediately 
before  0800  LT  and  continue  until  about  1800  LT.  Around  noontime, 
scintillations  attain  maximum  daytime  level  which  can  be  as  large 
as  13  dB  peak  to  peak  at  41  MHz  and  a few  dB  even  at  140  MHz. 

On  some  occasions,  the  nighttime  F-region  scintillation  continues 
through  the  early  morning  hours  to  merge  with  the  daytime  scintil- 
lation. On  these  days  it  is  impossible  to  determine  the  onset 
time  of  daytime  scintillation.  On  some  other  occasions,  the  onset 
of  post-sunset  F-region  scintillation  takes  place  before  the 
daytime  scintillation  event  has  ceased  completely. 

Figure  1 indicates  the  time  intervals  over  which  both  the  NOAA 
oblique  VIIF  radar  detected  electrojet  irregularities  and  scintil- 
lations were  recorded  at  Huancayo.  On  a few  days  scintillation  measure 
ments  at  41  MHz  only  were  available  and  on  some  others  scintil- 
lation measurements  at  140  MHz  only  could  be  made.  The  shading 
for  these  different  measurements  have  been  indicated  in  the 
diagrams.  The  breaks  in  the  bar  diagrams  indicate  either  no 
observations  were  made  during  the  period  or  no  satellite  trans- 
mission was  available.  During  the  period  February  1-8,  1975, 
scintillations  are  observed  only  when  50  MHz  radar  is  also  able 
to  detect  electrojet  irregularities.  In  view  of  the  higher 
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frequency,  scintillations  at  140  MHz  are  detected  slightly  after 
the  radar  detects  the  presence  of  electrojet  irregularities. 
Similarly,  the  cessation  of  scintillation  activity  at  140  MHz 
occurs  earlier  than  the  cessation  of  radar  echo.  We  shall 
later  provide  estimates  of  irregularity  amplitudes  required  to 
give  measurable  scintillation  at  140  MHz.  This  will  indicate 
that  the  irregularity  amplitudes  need  to  be  larger  than  some 
threshold  level  to  give  measurable  140  MHz  scintillation. 

It  is  interesting  to  note  that  on  February  4,  the  electrojet 
irregularities  are  absent  during  the  major  part  of  the  day  as 
indicated  by  the  radar  measurements  and  scintillations  are  also 
conspicuous  by  their  absence.  The  weak  irregularities  detected 
by  the  radar  between  1500  and  1800  LT  are  not  sufficiently  strong 
to  cause  140  MHz  scintillation.  Scintillation  at  140  MHz 
observed  on  February  1 immediately  prior  to  2000  LT  without  any 
associated  radar  echo  is  attributed  to  the  post-sunset  F region 
irregularities.  On  this  day,  the  nighttime  scintillation  at 
140  MHz  occurred  at  1945  LT  and  continued  until  0445  LT  on  the 
next  day.  Between  February  18  and  March  5,  1975,  scintillation 


measurements  could  be  made  at  the  lower  frequency  of  41  MHz  and 
it  may  be  noted  that  the  time  synchronism  between  the  radar  echo 
and  scintillation  is  excellent.  The  weak  irregularity  amplitudes 
which  become  detectable  with  the  radar  give  rise  to  immediate 
scintillation  effects  at  the  lower  frequency. 

The  comparative  study  of  radar  echoes  and  VHF  scintillation 
illustrated  by  Figure  1 indicate  that  daytime  scintillation  is 
caused  by  electrojct  irregularities  and  that  the  irregularities 
with  wavelengths  of  the  order  of  a km  (corresponding  to  Fresnel 
dimensions)  as  explored  by  scintillation  measurements  co-exist 
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with  the  short-wavelength  (3m)  irregularities  detected  by  the 
50  MHz  radar. 


We  shall  next  illustrate  the  behavior  of  q-type  sporadic-E 
(E  ) echoes  on  February  20  in  relation  to  41  MHz  scintillation 
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and  VII F radar  echo.  Figure  2 shows  that  both  the  scintillation 

activity  and  radar  echo  corresponding  to  type  2 ceased  at  1820 

LT.  The  bottom  panel  of  this  figure  shows  the  ionogram  sequence 

between  1 800  and  1830  LT.  The  presence  of  E echoes  in  the 

first  two  frames  may  be  observed  but  in  the  final  frame  obtained 

at  1830  LT,  Es^  is  seen  to  have  disappeared.  Such  a correlation 

between  these  events  is  found  to  exist  on  other  days  also.  It 

shows  that  VHF  scintillation,  type  2 radar  echo  and  E echoes 
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appear  and  disappear  together.  An  extensive  study  of  echoes 

in  relation  to  type  2 radar  echoes  has  been  recently  carried 
out.  Our  results  establish  the  fact  that  type  2 irregularities 
are  synchronized  not  only  with  the  appearance  of  irregularities 
at  30-80m  wavelengths  as  viewed  by  E echoes  on  ionograms  but 
with  kilometer  wavelength  irregularities  as  well  which  cause 
41  MHz  scintillations.  Thus  this  correlative  study  indicates 
that  type  2 irregularities  cover  about  a 3-decade  range  in 
wavelength . 

The  radar  observations  provided  information  on  the  nature 
of  electrojet  irregularities  observed  at  50  MHz.  It  has  been 
established  that  both  type  1 and  type  2 are  equally  important 
in  giving  50  MHz  radar  backscatter  echoes  whereas  at  frequencies 
higher  than  50  MHz,  type  1 irregularities  predominate  and  radar 
echoes  at  lower  frequencies  are  dominated  by  type  2 irregular- 
ities. The  present  50  MHz  radar  observations  yield  the  time 
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evolution  of  type  1 and  type  2 irregularities  and  also  provide 
the  strength  of  either  type  of  echo  simply  as  weak  or  strong. 
Figure  3 shows  the  time  variation  of  scintillation  index  SI  at 
140  MHz  and  the  evolution  of  type  1 and  type  2 irregularities 
as  observed  by  the  50  MHz  radar.  The  levels  1 and  2 marked  in 
NOAA  radar  plots  represent  respectively  the  presence  of  weak 
or  strong  type  2 irregularities  whereas  levels  3 and  4 correspond 
respectively  to  weak  and  strong  type  1 irregularities.  On 
February  2,  140  MHz  scintillation  becomes  discernible  at  0915  LT 
when  the  radar  echo  reveals  the  presence  of  strong  type  2 and 
weak  type  1 echo.  At  1210  LT  when  the  radar  echo  indicates  the 
cessation  of  type  1 irregularities  and  presence  of  type  2 
irregularities  only,  the  scintillation  index  does  not  show  any 
change.  This  implies  that  the  presence  or  absence  of  type  1 
irregularities  do  not  introduce  any  marked  variation  of  the 
scintillation  level  indicating  thereby  that  VHF  scintillation 
activity  is  related  to  type  2 irregularities. 

The  lower  data  set  obtained  on  February  4 shows  an  unusual 
day  with  no  electrojet  irregularities  and  consequently  no  scin- 
tillation during  the  major  part  of  the  day.  The  sporadic 
appearance  of  radar  echo  between  1600-1920  LT  corresponding 
mostly  to  weak  type  2 is  not  associated  with  scintillation 
activity  at  140  MHz.  The  absence  of  scintillation  activity  on 
that  day  as  well  as  any  type  2 irregularities  during  the  major 
part  of  the  day  may  be  related  to  a reversal  of  electric  field 
at  E-region  heights.  On  a normal  day,  the  E-region  electric 
field  near  the  magnetic  equator  is  eastward,  an  upward  Hall 
electric  field  is  generated  which  in  the  presence  of  an  upward 


plasma  density  gradient  generates  type  2 irregularities.  The 
generation  of  these  irregularities  is  inhibited  when  the  electric 
field  reverses  its  direction  to  become  westward.  The  reversal 


of  the  electric  field,  termed  as  count er -e 1 ectroj et , has  been 
diagnosed  from  a decrease  in  horizontal  magnetic  field  intensity 


and  has  been  shown  to  be  correlated  with  the  disappearance  of 


echoes  and  decrease  of  signal  intensity  in  transequator ial  VHF 


scatter  links  . 


During  the  development  phase  of  VHF  scintillations,  a scin- 
tillation index  (SI)  of  about  50%  at  41  MHz  typically  corresponds 
to  moderately  strong  type  2 echoes  at  meter  wavelengths.  The  SI 
index  of  50%  corresponds  to  the  theory  based  index  of  about 
0.25.  A knowledge  of  index  allows  us  to  compute  the  electron 
density  deviation  in  the  irregularities  of  realistic  assumptions 
can  be  made  regarding  the  irregularity  layer  thickness  (L),  layer 
height  (z)  and  irregularity  power  spectrum.  The  slant  range  (z) 
of  electrojet  region  to  Huancayo  for  ATS-6  observations  during 
the  period  was  about  110  km.  The  thickness  (L)  of  this  layer 
over  which  the  irregularities  are  assumed  to  be  distributed 
undiminished  in  magnitude  is  considered  as  5 km.  We  have  also 
assumed  that  the  3-dimensional  irregularity  power  spectrum  has 
a spectral  index  of  4 in  agreement  with  other  experimental  and 
theoretical  results.  The  outer  scale  of  the  irregularity  has 
been  taken  as  2 km  in  keeping  with  the  numerical  studies  of 
type  2 irregularities.  We  find  that  index  of  0.25  at  41  MHz 

may  be  observed  when  the  electron  density  deviation  of  the 
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irregularities  is  AN  = 6x10  m . This  corresponds  to  irreg- 
ularity amplitudes  of  6%  and  3%  respectively  for  amb i en t 
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ionization  density  of  10**  and  2xl0**m corresponding  to 
critical  frequencies  of  3 and  4 MHz.  At  140  MHz,  similar 
level  of  is  attained  only  if  the  irregularity  amplitudes 
increase  by  a factor  of  about  5.  Theoretical  studies  of  type  2 
irregularities  indicate  that  when  the  amplitude  of  type  2 
irregularities  in  the  tens  of  meter  wavelength  range  attains 
a level  of  about  4%  under  the  linear  theory  of  gradient  drift 
instability,  shorter  wavelength  irregularities  in  the  meterwave 
range  may  be  generated  through  non  linear  cascading  mechanism. 

Our  computations  of  irregularity  amplitudes  based  on  VHF  scin- 
tillation measurements  indicate  that  when  the  r.m.s.  irregularity 
amplitude  in  the  kilometer  wavelength  range  is  found  to  be  of 
the  order  of  5%,  the  50  MHz  radar  is  able  to  detect  simultaneously 
the  presence  of  type  2 irregularities.  The  highest  level  of 
index  at  41  MHz  during  daytime  is  found  to  be  of  the  order  of  0.5, 
which  is  found  to  correspond  to  an  electron  density  deviation 
AN  = 1.6xl0*°m~3  with  the  above  assumed  parameters.  This  value 

corresponds  to  r.m.s.  irregularity  amplitudes  of  16%  for  ambient 
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ionization  density  of  1C  m 

For  a power  law  type  of  irregularity  power  spectrum,  the 
scintillation  index  (S^)  is  expected  to  have  a power  law  variation 
with  frequency.  Assuming  that  such  a power  law  variation  exists^ 
we  defined  the  spectral  index  (n)  of  sc int i 1 1 at i on  as 

log  (S4  /S4  ) 

1 42 

n Log  (f2/f  ) 


where  S,  and  S.  are  the  scintillation  indices  (S.)  at  f,  and  f„ 
4 2 4 1* 
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respectively.  From  simultaneous  daytime  scintillation  measure- 
ments at  41  and  140  MHz  the  spectral  index  n is  obtained. 

Figure  4 indicates  the  variation  of  n with  S4  index  at  41  MHz. 

It  is  found  that  for  S^  varying  between  0.2  to  0.4,  n varies 
between  0.6  to  1.2.  For  weak  scattering  spectral  index  of 
scintillation  related  to  the  spectral  index  (p)  of  the  irreg- 
ularity power  spectrum  as  r)  = P-^-— . Thus  the  results  indicate 
that  during  the  development  or  decay  of  type  2 irregularities, 
the  irregularity  power  spectrum  is  very  flat  but  attains  a 
spectral  index  of  about  3 when  the  irregularities  are  fully 
developed.  It  should  be  emphasized  that  the  Fresnal  dimensions 
of  41  and  140  MHz  at  the  electrojet  height  corresponds  to 
0.5  km  and  1 km  respectively.  Thus  the  derived  spectral  indices 
of  the  irregularities  correspond  to  the  above  wavelength  range 
which  is  close  to  the  outer  scale  of  type  2 irregularities. 


3.  SUMMARY' 

The  study  reveals  that  daytime  VIIF  scintillations  are 


generated  when  Eso  echoes  appear  on  ionograms  and  50  MHz  radar 


detects  type  2 electrojet  irregularities.  Since  VHF  scintillations 


are  sensitive  to  kilometer  wavelength  irregularities,  E monitors 
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30-80m  wavelength  irregularities  and  50  MHz  radar  detects  3m  wavelength 

irregularities,  it  may  be  concluded  that  type  2 irregularities 
covering  the  wavelength  range  of  a few  meters  to  a few  kilometers 
co-exist.  We  have  also  found  that  during  periods  when  the  radar 


does  not  detect  any  electrojet  irregularity,  neither  E nor  VHF 


scintillation  is  observed. 


With  the  time  scale  of  our  measurements,  it  is  not  possible 


to  locate  the  time  period  when  non-linear  cascading  mechanism 
starts  to  generate  meter  wave  irregularities,  the  theoretically 
predicted  time  interval  being  of  the  order  of  a second. 

We  have  shown  that  when  type  2 meter  wave  irregularities 

and  VHF  scintillations  are  simultaneously  observed  the  irreg- 

ulatity  amplitudes  ranging  between  4%  to  16%  (for  ambient 
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ionization  density  of  10  m ) at  kilometer  wavelengths  are 
obtained.  Further  in  the  wavelength  range  of  0.5-1  km  corres- 
ponding to  the  Fresnal  dimensions  of  scintillation  measurements, 
the  spectral  index  of  the  irregularity  power  spectrum  indicates 
a variation  with  irregularity  amplitudes.  At  low  levels  of 
scintillation,  the  power  spectrum  is  very  flat  but  confirms  to 
a power  law  spectral  exponent  of  about  3 when  fully  developed. 
The  observed  evolution  of  irregularity  spectrum  corresponds  to 
wavelengths  in  the  vicinity  of  outer  scale  of  type  2 irreg- 
ularities which  is  typically  assumed  to  lie  in  the  range  of 
2-3  km. 
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SIGNAL  ANALYSIS  PACKAGE 

The  signal  analysis  package  described  in  a Final  Report 
AFGL- TR - 76 - 0 04  1 , has  been  modified  for  greater  ease  of  use. 
Following  is  a description  of  the  deck  setup  required  to 
operate  this  program  on  the  CDC  6600.  In  all  that  follows, 
optional  parameters  are  enclosed  in  square  brackets  [ ]. 


Deck  Setup 

Job  Card  Jobname ,CM1 10000 , Txxx,TPl . Job#  Name 
VSN,TAPEl=xxxxxx. 

REQUEST, TAPE1 , HY , S (xxxxxx/NORING) 
ATTACH, SIG, SI GNALYSISx3828, I D=C ANTOR . 
CALL, I=SIG,P=runoption. 

7/8/9 


Calibration  control  card[s] 
control  card[s] 

Title  card 

cont inuation 
run  s 


[Calibration  control  card[s]  ] 
[ cont  ro 1 card  [ s ] ] 

Title  card 


Run  option 

There  are  five  run  options: 


ONRED 
ON BLACK 
OF FRED 
OFF BLACK 
NOPLOT 


ON  or  OFF  for  on  or  offline  plot 
RED  or  BLACK  ink 


Calibration  Control  Cards 
Txx,Cxx,Fxx[Tj  , T£ , Rj , R2 ] [ , II]  [ , U ] 


Tape  number  xx  is  floating  point  (T37,  T45.1) 
Cal  number  (C3) 

File  number  (F2) 


r—mnmm m ■ k i - 
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T2  are  start  and  end  time  of  digitization  for  this  file 


Rj,  R2  are  start  and  end  time  (real  time  of  data  acquisition 
for  this  file) 


, R^  are  of  the  form  hhmmss 


Here  cals  follow 


Update  - update  cal  file 


Card  1 


Card  2-n  - 


Calibration  Cards  (H  option) 

Channel  number,  identification  (I2,A10)  - b2360MHZ 

digitized  count  level,  dB  level  free  format 
separated  by  blank  or  comma.  Need  not  be  one 
level  per  card.  Calibrations  are  entered  with 
high  signal  first  and  nondecreasing  signal 
strengths.  All  dB  levels  are  normalized  to 
0 dB  high  signal 

Exampl e : 

128,0,100,1,90,2  is  normalized  to 


Card  n+1 


90-2 

Calibrations  can  be  added  at  either  end, 

Examp  1 e : 

140,-2,135,-1 

128,0,100,1,90,2 

85,3  is  normalized  to 


i b 1 ankl 

l -9  J 


end  of 


f channel 
calibrations 


If  no  calibrations  are  entered 
(i.e.,  n ♦ 1 = 2 ) dB  = ievel. 
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Examples: 

1)  Existing  cal  file 
T3 7 , C2 , F 3 

2)  To  add  new  file  or  redefine  old  file  with  previously 
defined  calibrations 

T37,C2,F4=1800. 30,  1815 .05, 200, 400, U 

3)  To  use  new  calibrations  but  not  to  redefine 
T37.C2 , F3=1800, 1900, 200, 400, H 

cals 

It  should  be  noted  that  calibrations  are  catalogued  under 
tape  and  calibration  number.  Any  redefinition  of  cali- 
brations for  any  file  under  that  tape  and  calibration 
number,  redefines  calibrations  for  all  files. 

© 

Control  Cards 

Control  cards  are  free  form  concatenations  of  option 
selections  of  the  form 

Option  or  Option  = parameter  list 

(options  are  separated  by  commas) 

Parameter  lists  have  the  following  forms: 


a) 

OFF 

turn  the  option  off 

for 

all  channels 

b) 

none 

" " " on 

II 

II  II 

c) 

n 1 ' ‘ ' 

•’nk 

turn  the  option  off 

for 

channels  n^,  off 

for  other  channels 

d) 

h hmm . 

ss  time 

e) 

<dl 

number 

f) 

dB1# 

dt^  dB  levels  limits 

for 

dwell  analysis 

g) 

t1*  * * 

•^n 

times  for  dwell  analysis 


options 

parameter 
list  types 

meaning 

default 

END 

end  of  control  cards 

CHANS 

a,b,c 

restrict  analysis  to  1 

these  channels 

,2,. ..,10 

CDF 

a ,b , c 

CDF's  on  these  channels 

OFF 

DWELL 

a , b , c 

dwell  time  ahalysis 
for  these  channels 

OFF 

MSG 

a , b , c 

message  reliability 
statistics  for  these 
channel s 

OFF 

DATA 

a , b , c 

plot  data  for  these 
channe 1 s 

OFF 

SPECTRA 

a , b , c 

plot  smoothed  power 
spectra  for  these 
channels 

OFF 

RAW 

a , b , c 

plot  raw  power  spectra 
for  these  channels 

OFF 

AUTO 

a , b , c 

perform  autocorrelation 
for  these  channels 

OFF 

CROSS 

a,b,c* 

perform  cros s -corre 1 at i on 
for  these  channels 

OFF 

CORR 

e 

correlation  limit  for 
plotting  of  auto  and 
cross-correlations . 
CORR=0  - plot  all 
C0RR=2  - suppress  plots 

0.3 

SECONDS 

e 

number  of  seconds  (real 
time)  in  each  period  of 
analysis 

900 

PERIODS 

e 

number  of  periods  to 
analyze 

999 

REWIND 

rewind  data  tape 

SKIP 

b , e 

skip  d, (or  0)  files  on 
data  tape  after  skipping 
to  end  of  current  file 
if  not  already  there 

parameter 

options  list  types 


meaning 


default 


START 


CLEAR 


d 


real  time  start  of 
data  analysis 


current  real  time 
on  data  tape  or 
start  of  current 
data  file 


all  analysis  options 
and  all  channels  off 


NTH 


e use  every  n'th  point 


1 


DBS 


TIMES 


dB  range  for  dwell 
times  and  message 
reliability  analysis 

10  to 

+ 5 

time  ranges  for  dwell 
time  and  message 

*n 

tn+  1 

At 

reliability  analysis 

.25 

1 

.25 

1 

6 

. 5 

6 

10 

1 

10 

20 

2 

20 

40 

5 

40 

60 

10 

60 

t 

max 

all  times  in 

* Cross  correlation  channel  list  consists  of  sublists 
separated  by  semicolons.  The  first  channel  in  the 
sub  list  is  correlated  with  the  remaining  channels. 
The  last  sublist  may  be  terminated  by  a comma  or 
end  of  card. 

For  example: 

CROSS= 1 , 2 , 3 ; 2 , 3 , 4 
or 

C ROSS  = 1 ,2,3; 

2,3,4 

or 

C RO  S S = 1 ,2 
3 ; 2 , 3 , 4 

will  yield  cros s - co r re  1 at i ons  of  channel  1 vs  2 + 3 

and  2 vs  3 + 4 , while 

CROSS= 1 ,2,3 
2,3,4 

will  yield  cross- correlations  of  1 vs  2,3,  and  4. 
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TITLE  CARD 

A title  card  is  any  card  that  cannot  be  mistaken  for 
a calibration  card  or  control  card.  Care  should  be  taken 

| 

not  to  use  a control  card  option  word  as  the  first  word 
ofatitle. 

Processing  begins  after  each  title  card  is  recognized. 
Exampl e : 

■ i 

To  run  4 fifteen-minute  data  samples  on  channels 
1,3,4  on  data  tape  0S1 , file  three  starting  at  415. 
Statistics  wanted  are  only  print  statistics  including 
those  for  correlations.  Further,  one  hour  samples 
of  smoothed  spectra  and  autocorrelations,  decimated 
4 to  1 are  wanted  for  the  fourth  file  of  the  tape. 

Cross -correlations  are  wanted  only  for  channel  1 + 2 
versus  the  rest,  and  are  to  be  plotted  only  if 
greater  than  0.4.  There  are  5 channels  on  the  tape. 

The  calibrations  are  catalogued  as  tape  17,  calibration 
1,  files  2 and  3.  Plots  are  to  be  off  line,  red  ink. 

Deck 

CANTOR,  CM1 1000 ,TS00 . 

VSN , 7APE1=050001 . 

REQUEST, TAPE 1,HY,S. (OSO 00 1 /NOR I N6 ) 

ATTACH, SIG.SIGNALYSISX382 8 , I D=CANT0R . 

CALL, I=SIG,P=OFFRED. 

7/8/9 

T17.C1  ,F2 

PERICDS=4 .SKIP-2, START =415 , 

CHANS= 1 , 3 , 4 , CDF , DWELL , MSG 
AUTO, CROSS, CORRs 2, END 
TITLE  FOR  FIRST  SAMPLE  RUN 
T17.C1  ,F3 
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PERIODS =999 .SECONDS =3600 
NTH=4 .SKIP, CLEAR, SPECTRA, AUTO, 
CROSS=1,2,3,4,5;2,3,4,5,CORR=0.4,END 
TITLE  FOR  SECOND  SAMPLE  RUN 
6/7 / 8/9 


FOLLOWING  THIS  SECTION  ARE: 

CALIBRATION  ERROR  MESSAGES  (3  pages) 
ERROR  MESSAGES 

APPENDIX  1 - Catalogued  Procedures  on 
File  S I GNALY SISX3828 

APPENDIX  2 - Calibration  File  Format 

APPENDIX  3 - Control  Card  Decoder 

including 

Sample  Main  Program 

Functions 

Sample  Output 
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CALIBRATION  ERROR  MESSAGES  (2  of  3) 


’rV 

PoiL/uijy  ,,  . . ?'*W  h6? 

* iogihl*  JeProauction 


_ CEAim.IL 
36  C MHZ 


90 


— *ONREQ. 


OISPOSE.PLCT,' PL. 


LCaUNXZSAA  , IC=CANIQR,MR-1- 


COPYL.LGl ,LGO,  LG2. 

— UNLOAO,LGl,LGG. — 

ATTACH, PEN, ONI INEPEN. 

—HORARY,  PEN. 1 

ATTACH, TAPE  13 , CALI CX2  9h4, 10= CANTOR, PW  = CJC. 

— L&Z- 

EXTEND, TAPE10. 

— EX-I-T-. 

EXTEND, TAPE1C . 


_*QEJER£D.. 

DISPOSE, TAPE3C , * LR . 

attach, i r.i , rrf  nnTNy-’qAA  t t r-r. a ntit r , mr=  i . 

COPYL  , LG1  ,LGO, LG2. 

__iINLOAH.,L-GlJ,-LOC_. 

ATTACH, PEN, OFF  LINE  PEN . 

LTCRARY,  P.ETU 1 

ATT ACH,TAPF1Q , CALIGX2944, 10= CANTOR , PK=CJC . 

LG?. 

EX  TEND,  TAPE  10 . 

EX  I T . 

EXTEND, TAPE13 . 


-*OFFDLACK. ; 

DISPOSE , T APE3E  ,*LL. 

_A-T  T-ACH  ,L  G 1 , -Cr-CDO  IN  X 2-9  4 4 ^-I  &=C  A NT  0 R-,-MR=  L. 

COPYL , LG  1 , LGO, LG2 . 

UNLOAD, LG1 ,LGO. 

AT  T AC H, PEN, OFF  LI Nr PEN . 

LI ORARY , P EG . - 

ATT  ACH, TAPE  1C , CALI 0X2 944, I D= CANTOR, P9=CJC. 

EXTEND, TAPE  10 . 

-EXIT. 

EX  TEND, TAPE  13. 
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♦DNFU  flfX.  

DISPOSE, PLOT, * OL. 

attach. i r,  i . nn-.nRT>jy->q^  . tD-CA-NTOR.  . 

COPri_,LGl  ,LGO,LG2. 

_UNLOAD^J_Gl  ,-LGGL. , 

ATT ACH  ,PE N , ONL  INEPFN. 

I T R RASJC-4  PEtU 

ATTACH,T APE  13 , CALI PX2 94  4, I D=  CANTOR, PN=CJC. 

I r.P. 

EXTEND, TAPE  1C  . 

EX  IX* 

EXTEND, TAPF13  . 


— *NOPL-OX. 

AT T ACH, LG  1, DEC  DB INX 29 44 , ID=CANTOR, MR=1. 

_OOP_Y-U-LG-i-,.LGO  ,-LGP. 

UNLOAD, LG1, LGO.  ' 

-ATTACH, PEN,  ONI  INEPEN. : 

LIBRARY, PEN, 

— ATTACH  , T A PC  10  , CAL-I  PX2  944,  I D-C  ANT  OR,PW  = C JC  • 

LG2. 

— E-XXLN  0.,XA-EE  1 n , 

EXIT. 

-EXTEND, TAPE  10. 
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APPENDIX  2 
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Calibration  File  Format 

Catalogued  calibations  reside  on  a CDC  6600  mass  storage 
file.  See  the  CDC  FORTRAN  Extended  manual  for  general 
definitions  and  usage  details.  The  defining  parameters  for 
a particular  calibration  are  type  number  and  calibration  number, 
referred  to  hereafter  as  a calibration  id.  Each  calibration  id 
may  have  associated  with  it  up  to  twenty  files.  Calibration  id's 
are  numerically  indexed  in  the  main  file  index  in  the  order  of 
receipt,  excapt  that  additions  or  changes  to  an  existing 
calibration  id  overwrite  it.  Up  to  one  hundred  calibration  id's 
may  be  catalogued. 

The  following  is  a description  of  the  calibration  file 

Main  Index  Structure 

Pointer  to  control  record  1 
Pointer  to  calibration  subindex  1 


Pointer  to  control  record  n 
Pointer  to  calibration  index  n 
0 


0 
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s t rue  ture , 

Word 

1 


2n  - 1 
2n 

2n  + 1 


200 


AOS*  *90 


UNCLASSIFIED 


EMMANUEL  COLL  BOSTON  MASS  F/G  20/14 

STUDY  AND  ANALYSIS  OF  TOTAL  ELECTRON  CONTENT  AND  SCINTILLATION  — ETCCU) 
OCT  76  S BASU*  C J CANTOR*  J JOHANSON  F19628-76-C-0003 


SCIENTIFIC-1 


AFGL-TR-76-0260 


2*2 

^§>34490 


2-77 


I 


Control  Record  Format 


Word  1 2 

Tape#  Cal* 


3 7 

File  id  #1  (5  words) 


8 12 
File  id  #2 


Calibration  id  from  calibration  control  card 


Word  13 


97  98  102 

File  id  #20 


All  variables  in  the  control  record  are  strictly 
positive  floating  point  numbers  except  that  a 
null  file  id  consists  of  -999. 0's. 


File  id  Format 


Word  1 - File  # from  calibration  control  card 

2 - Start  time  (tape)  from  calibration  control  card 

3 - End  time  (tape)  from  calibration  control  card 

4 - Start  time  (real)  from  calibration  control  card 

5 - End  time  (real)  from  calibration  control  card 

(Time  formats  are  hhmm.ss) 

Calibration  Subindex  Structure 


Word 

1 Pointer  to  channel  1 calibration  record 

2 Pointer  to  channel  2 calibration  record 


10  Pointer  to  channel  10  calibration  record 

11  0 
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Calibration  Record  Format 


Word 


256 


dBs 

1 count 


dBs 

256  counts 


Word 


257 


258 


259  260 


0 if  signal  decreases 

with  counts. 

1 if  signal  increases 

with  count 
-1  see  Note  1 


#characters 
alpha  in 
259-260 


Alpha . from 
input 


d B ' s are  < for  all  counts  except  that  original 
input  calibrations  are  positive. 


Note  1.  A null  calibration  on  input  will  result  in 
setting  dB^=i.  Calibrations  should  not  be 
catalogued  this  way,  but  various  statistics 
on  uncalibrated  data  can  be  generated  this 
way. 


E xamp le : 

T999,C999,F999=t1,t2,r1,r2,H 

1137MHZ 

blank 

2256MHZ 

-9 

DATA .SPECTRA , START  = r ^ ,END 

QUICK  LOOK  AT  SPECTRA  OF  RAW  DATA. 


APPENDIX  3 

As  part  of  Signal  Analysis  Package,  a control  card 
decoder  was  written.  This  relies  heavily  on  two  functions 
NUM  and  IFIELD  which  are  of  more  general  interest. 

Taken  together,  these  functions  will  scan  a card,  stripping 
off  EDITOR  defined  sequence  numbers  if  desired,  and  return 
individual  fields  on  the  card  between  user  defined 
terminators . 

The  field  is  returned  either  as  alphabetic  array  or 
a floating  point  variable.  In  addition,  the  length  of  the 
field  and  the  terminator  is  available. 

Following  is  a sample  main  program,  the  functions, 
and  sample  output.  Usage  instructions  are  contained  in 
comments  to  the  functions. 
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